Prom second to first order transitions in a disordered quantum magnet 
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We study the spin-glass transition in a disordered quantum model. There is a region in the phase 
diagram where quantum effects are small and the phase transition is second order, as in the classical 
case. In another region, quantum fluctuations drive the transition first order. Across the first order 
line the susceptibility is discontinuous and shows hysteresis. Our findings reproduce qualitatively 
observations on LiHo a! Yi_ a ;F4. We also discuss a marginally stable spin-glass state and derive some 
results previously obtained from the real-time dynamics of the model coupled to a bath. 



The study of quantum effects on the properties of spin 
glasses is a subject of great experimental and theoret- 
ical interest. Spin-glass phases have been identified in 
systems such as mixed hydrogen-bonded ferroelectrics 
jy, the dipolar magnet LiHo 2: Yi_ :E F4 [|]|| or Sr-doped 
La2CuC>4 where quantum mechanics plays a funda- 
mental role. An important question is whether quantum 
spin glasses are qualitatively different from their classi- 
cal counterparts at low temperature. There is growing 
experimental evidence that the answer to this question 
is affirmative both in and out of equilibrium The 
thermodynamics of several models of disordered mag- 
netic systems has been investigated with various tech- 
niques. Mean-field-like models have been solved using 
the replica formalism in imaginary-time |p|-|l3|| and the 
Ising model in a transverse field has also been studied 
in finite dimensions Q-M. It is generally found that, 
in terms of a suitably defined quantum parameter T, a 
boundary T C (T) in the T — T plane separates spin-glass 
(SG) and paramagnetic (PM) phases. The transition line 
ends at a quantum critical point at T = 0, r c (0) above 
which the system is paramagnetic at all temperatures. 
In the case of the quantum spherical p-spin model, the 
real-time dynamics of the system coupled to a phonon 
bath was also investigated JlTj . In this boundary 
Fd(T) was found across which there is a dynamic phase 
transition from a PM state with equilibrium dynamics to 
a SG with non-stationary, aging, dynamics. 

In this paper we investigate in detail the equilibrium 
properties of this model. We find that a tricritical point 
(T*, r*) divides the line T C (T) in two parts. For T > T*, 
the SG transition is of second order and the behavior of 
the quantum system is qualitatively similar to that of the 
classical one. However, for T < T* quantum fluctuations 
drive the transition first order. The magnetic suscep- 
tibility is discontinuous and shows hysteresis across the 
first-order line. These findings reproduce qualitatively 
the observed behavior of LiHo x Yi„ x F4 in a transverse 
magnetic field The equations describing this sys- 

tem are non-linear and there is multiplicity of solutions in 
parts of the phase diagram. We found as a surprise that 



the usual criteria used to choose between them have to be 
reinterpreted in order to get physically meaningful solu- 
tions in the region T < T*. We also discuss the properties 
of solutions obtained through the use of the marginality 
condition, an approach recently applied to the study of 
quantum problems |K|Jlsl] . It is known from work on 
classical models Jl9| that the results of this approach are 
closely related to those obtained from the analysis of the 
real-time dynamics of the system. We explicitly show 
that this holds true in our quantum case. This unablcs 
us to identify the dynamical transition line Fd(T) and to 
derive certain properties of the non-equilibrium dynamics 
using the replica calculation. 

The Hamiltonian of the quantum p-spin spherical 
model is 

N N 

H[P, S ,J] = —J2Pi- E Jh-i„Sh-Si P , (1) 

i=l ii<...« p 

where Si is a scalar spin variable and the conjugated mo- 
menta Pi satisfy the commutation relations [Pi,Sj] = 
—ihSij. A Lagrange multiplier z enforces the spherical 

constraint l/^X)i=i( s ?) = 1- T ne interactions Ji 1 ...i p 
are chosen from a Gaussian distribution with zero mean 
and variance ]j = J 2 p\/(2NP~ 1 ). The model 

has glassy properties for all p > 2. The Hamiltonian 
([!]) may be interpreted in several ways. It represents a 
non-linear generalization of the quantum-rotor spin-glass 
models discussed in the literature [O]. It also describes 
a quantum particle moving in an N (eventually infinite) 
dimensional space in the presence of a random potential. 
Finally, its partition function is formally identical to that 
of a classical chain of "length" L = (3h embedded in an 
A-dimcnsional random environment. 

The equilibrium properties of the model are obtained 
using a replicated imaginary-time path integral formal- 
ism [b|. In the large N limit, the saddle-point evalu- 
ation of the partition-function allows us to define the 
order-parameter Q ab {r - t') = 1/N J2?=i ( r <(^KV)) 
where a, 6 = 1, . . . , n denote the replica indices and T 
the imaginary-time ordering operator. In terms of Q a b 
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the free-energy per spin reads 



F = lim — Y [Trln ( - n ((Mu 2 k + z) 
n-»o 2n y (3 ^ Y 1 " 1 



\ (3h 



n 
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ab 



drQl b (r)\ (2) 



where j3 — \/(kgT) is the inverse temperature, oj k = 
are the Matsubara frequencies, Q a b(^k) — 
Jq H dTQ ab {T)e luJkT and q d {u k ) = Q aa (u k ). ^.From here 
on we take J as the unit of energy, K/ J as the unit of time, 
and work with dimcnsionlcss quantities. Quantum fluc- 
tuations are controlled by the parameter V = h 2 /(JM). 
The classical limit of the model |2(J is recovered when 
r — > 0. The equilibrium solutions are determined by re- 
quiring that Qab(^k), parametrized according to different 
ansatze, be an extremum of F. In the following we con- 
centrate on the case p > 3. The phenomenology of the 
p = 2 case (l^]ll| is not as rich. 

For sufficiently high T and/or T, thermal and/or quan- 
tum fluctuations destroy the SG phase and the system is 
in the PM phase. The free-energy is then extremal for 
Qab(T~) — 5 ab qd{r). Its Fourier transform is the solution 
of the equation 



q~d(uk) = [^/r + z - s(wfe)] 



(3) 
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with S(wfe) = p/2 J Q drq d ~ {r)e l ^ kT and z is determined 
from qd(0) = 1. The above equation is non-linear and 
may have several solutions, some of which may be spuri- 
ous. We solved Eq. (||) numerically for p = 3. We found 
that for T > T* w 1/6 there is only one solution, irre- 
spective of the value of T. However, for T < T*, three 
solutions coexist in a finite region of the T — T plane 
(not including the T=0 axis). One of them is unstable 
and can be discarded from the start. We discuss below 
how to choose the physical solution between the remain- 
ing two. In the SG phase, inspired by the classical case 
p(| , we searched for one-step RSB solutions of the form 
Qab{r) = q' d {r)S ab + qEA^ab, where e ab — 1 if a and b 
belong to the same diagonal block of size mxm and zero 
otherwise, and q' d {r) = qd(j) — qea- The diagonal part, 
qd{i~), the breaking point, m, and the Edwards- Anderson 
order parameter, (?ea, are determined by extremizing F. 
We find 



q' d {"k) = H/T + z> (£'M - E'(0))] 1 , 
where E'(r) = p/2{qV-\ T ) - q^ 1 ), z' = P /2f3mq EA 



(4) 
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Xp)/x p and 



m = Tx 



p yj2/(p(l + x p ))q E Z 



p/2 



(5) 



The parameter x p , solution of an algebraic equation, de- 
pends on p only, and we found X3 = 1.817. The con- 
dition qd{0) = 1 now yields an equation for the break- 
ing point of the form m ee /i T (r). Solutions of Eq. (@) 



exist only for T < r max (T). Above this value, quan- 
tum fluctuations destroy the SG phase. We found that 
the function ht(F) has two real branches in the interval 
< r < r max (T). The physical values of m are on the 
lowest branch which verifies /ir(0) = Ti c i ass (T), the clas- 
sical breaking point parameter. The situations above and 
below T* are different. For T > T* (but lower than the 
classical transition temperature), m max = m(r max ) = 1, 
its largest possible value. For T < T* , instead, m max < 1. 
In both cases, <7ea is finite at r max (T). We discuss 
below the consequences of these facts. We found that 
limT^o ?7i max (T) = 0, implying that replica symmetry is 
restored at the quantum critical point as in the model dis- 
cussed in reference 10 1. All these conclusions, obtained 
from the numerical analysis of the p = 3 case, also follow 
from an approximate analytical solution of the equations 
for arbitrary p > 3 |21[ . 




FIG. 1. Free-energies of the different PM (solid lines) and 
SG (symbols) phases above (a) and below (b) T* . The inset 
in panel (b) shows in detail the crossing of the free-energies 
at the critical point for T <T* . 

PM solutions exist throughout the T — T plane. The 
free-energies of the different states must thus be com- 
pared in order to construct a phase diagram. Figs. 0(a) 
and (b) show the T-dependence of the PM and SG free- 
energies for the case p = 3 computed from Eq. (||) for two 
temperatures, above and below T* . Solid lines and sym- 
bols represent the PM and SG solutions, respectively. 
The curves end at the point where the corresponding 
solution disappears. It may be seen that for T > T* 
the free-energies of the two states intersect precisely at 
T C (T) = r max (T): the SG solution does not extend be- 
yond the transition point and no hysteresis is expected. 
Below the critical point, F§g > Fpm meaning that the 
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SG solution maximizes the free-energy. This is the usual 
situation encountered in replica theories of classical spin 
glasses. As in the classical case, qea is discontinuous at 
the transition. The latter is nevertheless of second order 
because m = 1 at T c and, therefore, the effective num- 
ber of degrees of freedom participating in the transition 
(1 — m)(/EA —> at T c . There is no latent heat and the 
linear susceptibility is continuous. 
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FIG. 2. Static (thin lines) and dynamic (thick lines) phase 
diagrams of the p-spin model for p = 3. Solid and dashed lines 
represent second and first order transitions, respectively. 

The situation is more involved below T*. To start with, 
one has to choose between the two PM solutions labeled 
PMi and PM 2 in Fig.g. Naively, one would choose the 
solution with the lowest free-energy, i.e., PM2. How- 
ever, this solution has unphysical properties. As shown 
in Fig. 0(b), its free-energy never intersects that of the 
SG phase. Both the free-energy and the susceptibility 
diverge as T — > 0. Furthermore, this solution disappears 
at a finite value of T (not shown in the figure) and cannot 
thus be reached starting from Y — 00. On the other hand, 
it is clear that the ground-state of Hamiltonian ([!]) must 
have finite susceptibility and energy. We thus conclude 
that PM 2 is a spurious solution and that PM4 has to 
be chosen even if its free-energy is higher |E| . The free- 
energies of the SG and PMi states cross at T c < r max as 
shown in the inset in Fig. 0(b). In the low temperature 
phase, Fsg < Fpm, the opposite of what we found for 
T > T*. The SG and PM solutions extend beyond the 
point where they cross. There is a region of phase co- 
existence and hysteresis effects are thus expected in the 
behavior of observables. 

Since now c/ea and m are discontinuous at T c the ther- 
modynamic transition is first order with latent heat and 
discontinuous susceptibility (see below). The phase dia- 
gram resulting from this analysis is represented in Fig.^ 
(thin lines). The fiat section is the first-order line. We 
have computed the Edwards-Anderson order parameter 
and the susceptibility, x — Jo ' dT[qd(r) — (1 — to^ea], 
as functions of Y for the p = 3 model. The results are 
displayed in Fig.||. The susceptibility has a cusp at T c 



for T > T* and a discontinuity for T < T*. The dotted 
lines correspond to the regions of metastability. Their 
end points give the amplitude of the ideal hysteresis cy- 
cle. The importance of quantum fluctuations may be 
appreciated from the fact that half way from the transi- 
tion the order parameter is already reduced by a factor of 
two. It can be shown analytically that the spectrum 
of magnetic excitations at T = is gaped both in the PM 
and SG phases for all Y ^ (see below, however). Conse- 
quently, the latent heat vanishes exponentially as T — > 0. 
Since it also vanishes at T* , it must have a maximum at 
some intermediate temperature. 

First order quantum transitions were also found in two 
other models, the fermionic SK-like spin-glass model || 
and a p-spin model in a transverse field pj]] . In contrast, 
the SG transitions of the Heisenberg EA model and of the 
SK model in a transverse field are known to be second 
order |J. This is also true in finite dimensions [fj^-[l6||. 
Early experiments on LiHo a; Yi_ a; F4 gave some indica- 
tions that the second order SG transition seen above 25 
mK, might become first order at lower temperatures [0. 
More recently, hysteresis effects have been observed in 
this system as a function of the transverse field 0, giv- 
ing further support to this idea. The model that we study 
here captures this phenomenology. 
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FIG. 3. Magnetic susceptibility (a) and Edwards- Anderson 
order parameter (b) of the p=3 model. 

We discuss next the consequences of the use of the 
marginality condition jlitl rather than thermodynamics 
in the determination of the breaking point. In this ap- 
proach it is not required that F be an extremum with 
respect to m but that the SG phase be marginally sta- 
ble. This implies that the "replicon" eigenvalue A must 
vanish throughout the low-temperature phase. The cal- 
culation of A (2l) is analogous to the classical one [po[ . 
The result is 
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A = 



q d (0) +/3<2ea(to- 1) 



92(0) 

/3 2 



/? 2 <z| a (to - 1) + q d (0)PQEA(m - 2) 



P-2 



^ PO - 1)?EA 

The value of m follows from the equation A = 0. 
combined with the equation <5F/<5<7ea = 0, yields 



p 2 

(6) 
that 

(7) 



m = T(p-2)y/2/{p(p- 

Notice that this expression is equivalent to Eq. (j^) 
with the substitution x p — * (p — 2). Interestingly 
enough, Eq. (^) is identical to the equation found for the 
fluctuation-dissipation theorem (FDT) violation parame- 
ter, X, in the real-time dynamical calculation J!?]]. More- 
over, the static and dynamical equations for gEA are also 
identical which implies that m — X . The coincidence be- 
tween the values of X and m for the marginal SG state 
has been noticed several times for classical models. This 
is the first explicit evidence of its validity in a quantum 
problem. X is related to the effective temperature |2j| 
of the system, T e g = X~ 1 T, where T is the tempera- 
ture of the thermal bath it is in contact with. Values of 
X ^= 1 signal the presence of a non-stationary dynamics 
and of FDT violations. The fact that f3X — (3m — > const 
when T — > shows that a non-trivial T e g is generated 
even when the temperature of the bath vanishes. We 
have also shown analytically |2lJ that the internal energy, 
computed from U = d(fJF)/d(3 at constant m, coincides 
with the long-time limit of the energy per spin as ob- 
tained from dynamics The dynamic transition line 
I^i (T) may be thus identified as the boundary of the re- 
gion in the T—T plane where the marginally stable SG ex- 
ists. Below this line, the dynamics of the system becomes 
non-stationary and FDT violations set in. The dynamic 
phase diagram for p = 3 is shown in Fig.|^ (thick lines). 
As in the equilibrium case, m is discontinuous across the 
dashed line. Td lies always above T c suggesting that the 
equilibrium state can never be reached dynamically start- 
ing from an initial state in the PM phase. The two lines 
are extremely close to each other for T ~ T* . Within the 
accuracy of our calculations we cannot assert whether 
they precisely touch at T*, an intriguing possibility. In 
the region T < T* , m varies continuously along T^(T) 
and vanishes at the quantum critical point. This has a 
consequence of potential interest for experiment: FDT 
violations are predicted to appear suddenly rather than 
gradually as is crossed coming from the high T region 
for T < T*. The stationary part of the time-dependent 
susceptibility in the SG phase can be calculated by an- 
alytic continuation of q' d {ujk). It may be shown that the 
excitation spectrum of the marginal SG state is gapless 
pl| . Furthermore, x"(w) may be calculated exactly for 
T, uj -> 0. The result is 
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Pip ~ 1) 
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A linear excitation spectrum has also been found in the 
case of the Heisenberg spin glass model JTo|] . However, 
a gapless spectrum is not a consequence of Goldstone's 
theorem here as our model does not have any continuous 
symmetry. A more extended discussion of our results will 
be presented in a forthcoming paper |Q . 
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